Inflation: Homogeneous Limit 



V. Mukhanov 
ASC, Department of Physics, LMU, Munich 

February 5, 2008 



Abstract 

I review the motivation for the early-time cosmic acceleration stage in 
expanding universe and discuss simple inflationary scenarios. Preheating 
and reheating are considered in great detail. This is a sample chapter from 
my book "Physical Foundations of Cosmology" published by Cambridge 
University Press (2005). 

Matter is distributed very homogeneously and isotropically on scales larger 
than a few hundred Mpc. The microwave background gives us a "photograph" 
of the early universe, which shows that at recombination the universe was ex- 
tremely homogeneous and isotropic (with accuracy ~ 10~ 4 ) on all scales up to 
the present horizon. Given that the universe evolves according to the Hubble 
law, it is natural to ask which initial conditions could lead to such homogeneity 
and isotropy. 

To obtain an exhaustive answer to this question we have to know the exact 
physical laws which govern the evolution of the very early universe. However, 
as long as we are interested only in the general features of the initial condi- 
tions it suffices to know a few simple properties of these laws. We will assume 
that inhomogeneity cannot be dissolved by expansion. This natural surmise 
is supported by General Relativity (see Part II of this book for details). We 
will also assume that nonperturbative quantum gravity does not play an es- 
sential role at sub-Planckian curvatures. On the other hand, we are nearly 
certain that nonperturbative quantum gravity effects become very important 
when the curvature reaches Planckian values and the notion of classical space- 
time breaks down. Therefore we address the initial conditions at the Planckian 
time U=tpi ~ 1CT 43 s. 

In this chapter we discuss the initial conditions problem we face in a de- 
celerating universe and show how this problem can be solved if the universe 
undergoes a stage of the accelerated expansion known as inflation. 

1 Problem of initial conditions 

There are two independent sets of initial conditions characterizing matter: a) 
its spatial distribution, described by the energy density e (x) and b) the initial 
field of velocities. Let us determine them given the current state of the universe. 



Homogeneity, isotropy (horizon) problem. The present homogeneous, 
isotropic domain of the universe is at least as large as the present horizon scale, 
cto ~ 10 28 cm. Initially the size of this domain was smaller by the ratio of 
the corresponding scale factors, aj/ao- Assuming that inhomogeneity cannot be 
dissolved by expansion, we may safely conclude that the size of the homogeneous, 
isotropic region from which our universe originated at t = tj was larger than 

k ~ ct — . (1) 
It is natural to compare this scale to the size of a causal region l c ~ cti : 

lc ti a 

To obtain a rough estimate of this ratio we note that if the primordial radiation 
dominates at tj ~ tpi, then its temperature is Tpi ~ 10 32 K. Hence 

(ai/ao) ~ (To/Tpi) ~ 1(T 32 

and we obtain 

Thus, at the initial Planckian time, the size of our universe exceeded the causal- 
ity scale by 28 orders of magnitude. This means that in 10 84 causally discon- 
nected regions the energy density was smoothly distributed with fractional vari- 
ation not exceeding Se/s ~ 10~ 4 . Because no signals can propagate faster than 
light, no causal physical processes can be responsible for such an unnaturally 
fine-tuned matter distribution. 

Assuming that the scale factor grows as some power of time, we can use an 
estimate a/t ~ a and rewrite formula J2j as 

— ~ — . (4) 

lc ao 

Thus, the size of our universe was initially larger than that of a causal patch 
by the ratio of the corresponding expansion rates. Assuming that gravity was 
always attractive and hence was decelerating the expansion, we conclude from 
iJH that the homogeneity scale was always larger than the scale of causality. 
Therefore, the homogeneity problem is also sometimes called the horizon prob- 
lem. 

Initial velocities (flatness) problem. Let us suppose for a minute that 
someone has managed to distribute matter in the required way. The next ques- 
tion concerns initial velocities. Only after they are specified is the Cauchy 
problem completely posed and can the equations of motion be used to predict 
the future of the universe unambiguously. The initial velocities must obey the 
Hubble law because otherwise the initial homogeneity is very quickly spoiled. 
That this has to occur in so many causally disconnected regions further compli- 
cates the horizon problem. Assuming that it has, nevertheless, been achieved, 
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we can ask how accurately the initial Hubble velocities have to be chosen for a 
given matter distribution. 

Let us consider a large spherically symmetric cloud of matter and compare 
its total energy with the kinetic energy due to Hubble expansion, E k . The total 
energy is the sum of the positive kinetic energy and the negative potential energy 
of the gravitational self- interaction, E p . It is conserved: 

E tot = E k + E P = E k + E P^ 

Because the kinetic energy is proportional to the velocity squared, 

E k = El (d t /d ) 2 

and we have 

E\ E\ E k \a J ■ 

Since E% ~ and a /di < 10" 28 , we find 

rptot 

=|r < "r 66 . (6) 

This means that for a given energy density distribution the initial Hubble ve- 
locities must be adjusted so that the huge negative gravitational energy of the 
matter is compensated by a huge positive kinetic energy to an unprecedented 
accuracy of 10 _54 %. An error in the initial velocities exceeding 10 _54 % has a 
dramatic consequence: the universe either recollapses or becomes "empty" too 
early. To stress the unnaturalness of this requirement one speaks of the initial 
velocities problem. 

Problem 1. How can the above consideration be made rigorous using the 
Birkhoff theorem? 



In general relativity the problem described can be reformulated in terms of 
the cosmological parameter ft (t). Using the definition of fl (t) we can rewrite 
Friedmann equation as 

(Ha) 



and hence 



cu - 1 = (no - 1) = (n - 1) ( < io- 56 . (8) 



(Ha) 



O; 



Note that this relation immediately follows from JSJ if we take into account that 
n = \E P \ /E k . We infer from (JSJ that the cosmological parameter must initially 
be extremely close to unity, corresponding to a flat universe. Therefore the 
problem of initial velocities is also called the flatness problem. 
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Initial perturbation problem. One further problem we mention here for 
completeness is the origin of the primordial inhomogeneitics needed to explain 
the large-scale structure of the universe. They must be initially of order Se/e ~ 
10 -5 on galactic scales. This further aggravates the very difficult problem of 
homogeneity and isotropy, making it completely intractable. We will see later 
that the problem of initial perturbations has the same roots as the horizon 
and flatness problems and that it can also be successfully solved in inflationary 
cosmology. However, for the moment we put it aside and proceed with the 
"more easy" problems. 

The above considerations clearly show that the initial conditions which led 
to the observed universe are very unnatural and non-generic. Of course, one can 
make the objection that naturalness is a question of taste and even claim that 
the most simple and symmetric initial conditions are "more physical." In the 
absence of a quantitative measure of "naturalness" for a set of initial conditions 
it is very difficult to argue with this attitude. On the other hand it is hard 
to imagine any measure which selects the special and degenerate conditions 
in preference to the generic ones. In the particular case under consideration 
the generic conditions would mean that the initial distribution of the matter is 
strongly inhomogeneous with 5s /e > 1 everywhere or, at least, in the causally 
disconnected regions. 

The universe is unique and we do not have the opportunity to repeat the 
"experiment of creation" many times. Therefore cosmological theory can claim 
to be a successful physical theory only if it can explain the state of the ob- 
served universe using simple physical ideas and starting with the most generic 
initial conditions. Otherwise it would simply amount to "cosmic archaeology," 
where "cosmic history" is written on the basis of a limited number of hot big 
bang remnants. If we are pretentious enough to answer the question raised by 
Einstein, "What really interests me is whether God had any choice when he 
created the World," we must be able to explain how a particular universe can 
be created starting with generic initial conditions. The inflationary paradigm 
seems to be a step in the right direction and it strongly restricts "God's choice." 
Moreover, it makes important predictions which can be verified experimentally 
(observationally) , thus giving cosmology the status of a physical theory. 

2 Inflation: main idea 

We have seen so far that the same ratio, ch/clq, enters both sets of independent 
initial conditions. The large value of this ratio determines the number of causally 
disconnected regions and defines the necessary accuracy of the initial velocities. 
If gravity was always attractive, then hi /ho is necessarily larger than unity 
because gravity decelerates an expansion. Therefore, the conclusion cn/ao 1 
can be avoided only if we assume that during some period of expansion gravity 
acted as a "repulsive" force, thus accelerating the expansion. In this case we can 
have hi/ ho < 1 and the creation of our type of universe from a single causally 
connected domain may become possible. A period of accelerated expansion is a 
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a 



decelerated Friedmann expansion 




Figure 1: 



necessary condition, but whether is it also sufficient depends on the particular 
model where this condition is realized. With these remarks in mind we arrive 
at the following general definition of inflation: 

Inflation is a stage of accelerated expansion of the universe when gravity acts 
as a repulsive force. 

Figure 1 shows how the old picture of a decelerated Friedmann universe is 
modified by inserting a stage of cosmic acceleration. It is obvious that if we do 
not want to spoil the successful predictions of the standard Friedmann model, 
such as nucleosynthesis, inflation should begin and end sufficiently early. We 
will see later that the requirement of the generation of primordial fluctuations 
further restricts the energy scale of inflation; namely, in the simple models 
inflation should be over at i/ <~ 10~ 34 to 10~ 36 s. Successful inflation must also 
possess a smooth graceful exit into the decelerated Friedmann stage because 
otherwise the homogeneity of the universe would be destroyed. 

Inflation explains the origin of the Big Bang; since it accelerates the ex- 
pansion, small initial velocities within a causally connected patch become very 
large. Furthermore, inflation can produce the whole observable universe from 
a small homogeneous domain even if the universe was strongly inhomogeneous 
outside of this domain. The reason is that in an accelerating universe there 
always exists an event horizon. It has size 

^max &max 

r e (t) = a(t) ( - = a{t) ^ (9) 
J a j aa 

t a(t) 

The integral converges even if a max — * oo because the expansion rate d grows 
with a. The existence of an event horizon means that anything at time t a dis- 
tance larger than r e (t) from an observer cannot influence that observer's future. 
Hence the future evolution of the region inside a ball of radius r e (t) is completely 
independent of the conditions outside a ball of radius 2r e (t) centered at the same 
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place. Let us assume that at t = ti matter was distributed homogeneously and 
isotropically only inside a ball of radius 2r e (ti) (Fig. 2). Then an inhomogene- 
ity propagating from outside this ball can spoil the homogeneity only in the 
region which was initially between the spheres of radii r e (ti) and 2r e (ti) . The 
region originating from the sphere of radius r e (ti) remains homogeneous. This 
internal domain can be influenced only by events which happened at U between 
the two spheres, where the matter was initially distributed homogeneously and 
isotropically. 

The physical size of the homogeneous internal region increases and is equal 

to 

r h (t f )=r e (ti)^- (10) 

CLi 

at the end of inflation. It is natural to compare this scale with the particle 
horizon size, which in an accelerated universe can be estimated as 

t a 

, n f dt , . f da a(t) . . , . 

r p (t)=a(t) -=a(t) ^r e (U) , 11 

J a J aa a t 

ti o.i 

since the main contribution to the integral comes from a ~ . At the end of in- 
flation r p (tf) <~ rh(tf) , that is, the size of the homogeneous region, originating 
from a causal domain, is of order the particle horizon scale. 

Thus, instead of considering a homogeneous universe in many causally dis- 
connected regions, we can begin with a small homogeneous causal domain which 
inflation blows up to a very large size, preserving the homogeneity irrespective 
of the conditions outside this domain. 

Problem 2. Why does the above consideration fail in a decelerating universe? 

The next question is whether we can relax the restriction of homogeneity 
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on the initial conditions. Namely, if we begin with a strongly inhomogeneous 
causal domain, can inflation still produce a large homogeneous universe? 

The answer to this question is positive. Let us assume that the initial energy 
density inhomogeneity is of order unity on scales ~ -H,- j that is, 

'A .1^^1.0(1), (12) 

£ ) t . £ Oi £ Chi 

where V is the spatial derivative with respect to the comoving coordinates. At 
t 3> ti, the contribution of this inhomogeneity to the variation of the energy 
density within the Hubble scale H (t)~ can be estimated as 

£ ) t £ a (tj a(t) 

where we have assumed that |Ve| je does not change substantially during ex- 
pansion. This assumption is supported by the analysis of the behavior of linear 
perturbations on scales larger than the curvature scale ti~ x (see chapters 7 and 
8). It follows from Ijl3|l that if the universe undergoes a stage of acceleration, 
that is, a (t) > d* for t > ti, then the contribution of a large initial inhomogene- 
ity to the energy variation on the curvature scale disappears. A patch of size 
H^ 1 becomes more and more homogeneous because the initial inhomogeneity is 
"kicked out" : the physical size of the perturbation, oc a, grows faster than the 
curvature scale, H^ 1 = a /a, while the perturbation amplitude does not change 
substantially. Since inhomogeneities are "devalued" within the curvature scale, 
the name "inflation" fairly captures the physical effect of accelerated expansion. 
The consideration above is far from rigorous. However, it gives the flavor of the 
"no-hair" theorem for an inflationary stage. 

To sum up, inflation demolishes large initial inhomogeneities and produces a 
homogeneous, isotropic domain. It follows from l|13(l that if we want to avoid the 
situation of a large initial perturbation reentering the present horizon, ~ Hq 1 , 
and inducing a large inhomogeneity, we have to assume that the initial expansion 
rate was much smaller than the rate of expansion today, that is, dj /gsq <C 1. 
More precisely, the CMB observations require that the variation of the energy 
density on the present horizon scale does not exceed 10 -5 . The traces of an 
initial large inhomogeneity will be sufficiently diluted only if d^/do < 10~ 5 . 
Rewriting relation (JSJl as 

n = i + (n,-i)Q^ , (i4) 



we see that if |f2j — 1| ~ O (1) then 

n Q = i (15) 

to very high accuracy. This important robust prediction of inflation has a kine- 
matical origin and it states that the total energy density of all components of 
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matter, irrespective of their origin, must be equal to the critical energy density 
today. We will see later that amplified quantum fluctuations lead to the tiny 
corrections to CIq = 1 , which are of order 1CP 5 . It is worth noting that, in 
contrast to a decelerating universe where f2 (t) — ► 1 as t — > 0, in an accelerating 
universe (t) — > 1 as t — > oo, that is, O = 1 is its future attractor. 

Problem 3. Why does the consideration above fail for f2j = 0? 



3 How can gravity become "repulsive"? 



To answer this question we recall the Friedmann equation: 

a = — — G(e + 3p)a. 



(16) 



Obviously, if the strong energy dominance condition, e+3p > 0, is satisfied, then 
a < and gravity decelerates the expansion. The universe can undergo a stage 
of accelerated expansion with a > only if this condition is violated, that is, if 
e + 3p < 0. One particular example of "matter" with broken energy dominance 
condition is a positive cosmological constant, for which py = —£y and e + 
3p = —2sy < 0. In this case the solution of Einstein's equations is a de Sitter 
universe discussed in previous Chapters. For t ^> H^ 1 , the de Sitter universe 
expands exponentially quickly, a cx exp (H\t) , and the rate of expansion grows 
as the scale factor. The exact de Sitter solution fails to satisfy all necessary 
conditions for successful inflation: namely, it does not possess a smooth graceful 
exit into the Friedmann stage. Therefore, in realistic inflationary models, it can 
be utilized only as a zero order approximation. To have a graceful exit from 
inflation we must allow the Hubble parameter to vary in time. 

Let us now determine the general conditions which must be satisfied in a 
successful inflationary model. Because 



a 
a 



H 2 +H, 



(17) 



and a should become negative during a graceful exit, the derivative of the Hubble 

/H 2 grows toward the 



H 



constant, H, must obviously be negative. The ratio H 

end of inflation and the graceful exit takes place when 

H 2 . Assuming that H 2 changes faster than H , that is, 
the following generic estimate for the duration of inflation: 

tf ~ Hi/ H. 



H 



becomes of order 
< 2HH, we obtain 

(18) 



where Hi and Hi refer to the beginning of inflation. At t ~ tf the expression 
on the right hand side in equation (|17fl changes sign and the universe begins to 
decelerate. 
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Inflation should last long enough to stretch a small domain to the scale of 
the observable universe. Rewriting the condition d^/do < 10~ 5 as 



a% af 
hf d 



a,; Hi at 



a,f Hf d 



< 10" 



and taking into account that d//do should be larger than 10 28 , we conclude that 
inflation is successful only if 

£/ > 10 33 ^i 



Let us assume that 



H, 



H f 



<C Hf and neglect the change of the Hubble parameter. 



Then the ratio of the scale factors can be roughly estimated as 
af/ai ~ exp (Hitf) ~ exp {h 2 / H t ) > 10 33 . 



(19) 



Hence inflation can solve the initial conditions problem only if tf > 75 H ~ , that 
is, it lasts longer than seventy five Hubble times (e-folds). Rewritten in terms 
of the initial values of the Hubble parameter and its derivative, this condition 
takes the form 

i'; 

(20) 



Hf 



1 

< 75' 



Using the Friedmann equations, we can reformulate it in terms of the bounds 
on the initial equation of state 



< 10" 



(21) 



Thus, at the beginning of inflation the deviation from the vacuum equation of 
state must not exceed one percent. Therefore an exact de Sitter solution is a 
very good approximation for the initial stage of inflation. Inflation ends when 
s +p — e. 

Problem 4. Consider an exceptional case where H decays at the same rate 

as H 2 , that is, H = —pH 2 , where p = const. Show that for p < 1 we have 
power-law inflation. This inflation has no natural graceful exit and in this 
sense is similar to a pure de Sitter universe. 



4 How to realize the equation of state p « — s 

Thus far we have used the language of ideal hydrodynamics, which is an ade- 
quate phenomenological description of matter on large scales. Now we discuss a 
simple field-theoretic model where the required equation of state can be realized. 
The natural candidate to drive inflation is a scalar field. The name given to such 
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a field is the "inflaton." We saw that the energy-momentum tensor for a scalar 
field can be rewritten in a form which mimics an ideal fluid. The homogeneous 
classical field (scalar condensate) is then characterized by energy density 

e=\tf + V{<p), (22) 

and pressure 

P=~¥> 2 -V». (23) 

We have neglected spatial derivatives here because they become negligible soon 
after the beginning of inflation due to the "no-hair" theorem. 

Problem 5. Consider a massive scalar field with potential V — \m}p 2 , 
where m <C mpi, and determine the bound on the allowed inhomogencity 
imposed by the requirement that the energy density must not exceed the 
Planckian value. Why does the contribution of the spatial gradients to 
the EMT decay more quickly than the contribution of the mass term? 

It follows from equations (|22H and 123|) that the scalar field has the desired 
equation of state only if p 2 -C V (<p) . Because p = — e + p 2 , the deviation of 
the equation of state from that for the vacuum is entirely characterized by the 
kinetic energy, ip 2 , which must be much smaller than the potential energy V (ip) . 
Successful realization of inflation thus requires keeping p 2 small compared to 
V (ip) during a sufficiently long time interval, or more precisely, for at least 75 
e-folds. In turn this depends on the shape of the potential V (<p). To determine 
which potentials can provide us with inflation, we have to study the behavior 
of a homogeneous classical scalar field in an expanding universe. The equation 
for this field can be derived either directly from the Klein-Gordon equation or 
by substituting (|22|l and (|23[) into the conservation law e — ~3H (e +p). The 
result is 

Cp + 3H<p + V tV =0, (24) 

where V v = dV/dp. This equation has to be supplemented by the Friedmann 
equation: 

ff2= f G^+^O' (25) 

where we have set G — 1 and k — 0. We first find the solutions of equations 
H24|l and (|25|l for a free massive scalar field and then study the behavior of the 
scalar field in the case of a general potential V (p) . 

4.1 Simple example: V = -m 2 ip 2 . 

Substituting H from (|25|l into (|24|l . we obtain the closed form equation for ip, 
yi + Vl2^ (tp 2 +m 2 p 2 ) 1/2 p + m 2 p = 0. (26) 
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(27) 



This is a nonlinear second order differential equation with no explicit time de- 
pendence. Therefore it can be reduced to a first order differential equation for 
ip (ip) . Taking into account that 

. dtp 

equation (|26|l becomes 

dtp V127T (tp 2 + m 2 ip 2 ) ip + m 2 ip 
dtp tp 

which can be studied using the phase diagram method. The behavior of the 
solutions in tp — tp plane is shown in Figure 3. The important feature of this 
diagram is the existence of an attractor solution to which all other solutions con- 
verge in time. One can distinguish different regions corresponding to different 
effective equations of state. Let us consider them in more detail. We restrict 
ourselves to the lower right quadrant (tp > 0, ip < 0); solutions in the other 
quadrants can easily be derived simply by taking into account the symmetry of 
the diagram. 

Ultra-hard equation of state. First we study the region where \ip\ 3> mip. 
It describes the situation when the potential energy is small compared to the 
kinetic energy, so that tp 2 ^S> V. It follows from equations (1221) and 1|23|) that in 
this case the equation of state is ultra-hard, p ~ +e. Neglecting rrup compared 
to ip in equation (|27p. we obtain 

d'-P /To - • 



, Yl-Kib. (28) 
dp 



The solution of this equation is 



C exp (y 



(29) 

where C < is a constant of integration. In turn, solving l|29|l for tp (t) gives 

ip = const -= Int. (30) 

V 12-7T 
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Substituting this result into equation (|25[) and neglecting the potential term, we 
obtain 

~-h < 31 > 

It immediately follows that a oc t 1 / 3 and e cx a~ 6 in agreement with the ultra- 
hard equation of state. Note that the solution obtained is exact for a massless 
scalar field. According to (|29|) the derivative of the scalar field decays exponen- 
tially more quickly than the value of the scalar field itself. Therefore, the large 
initial value of \p>\ is damped within a short time interval before the field tp itself 
has changed significantly. The trajectory which begins at large \ip\ goes up very 
sharply and meets the attractor. This substantially enlarges the set of initial 
conditions which lead to an inflationary stage. 

Inflationary solution. If a trajectory joins the attractor where it is flat, at 
\ip\ ^> 1, then afterwards the solution describes a stage of accelerated expansion 
(recall that we work in Planckian units). To determine the attractor solution 
we assume that dtp /dtp ~ along its trajectory. It follows from 127fl that 

in 

(p at r ~ 7==, (32) 



and therefore 



TTl 

^atr (t) ^ (fit 7r==(^ - U) — r— (tf ~ t) , (33) 



where ti is the time when the trajectory joins the attractor and tf is the moment 
when ip formally vanishes. In reality, solution l|33|l fails well before the field p 
vanishes. 

Problem 6. Calculate the corrections to the approximate attractor solution 
(|32|) and show that 

p atI = — JL= (1 - ~ (ViW) 2 + o ( (ViW) 3 )) ■ (34) 



The corrections to (|32f) become of order the leading term when ip ~ O (1) , 
that is, when the scalar field value drops to the Planckian value or, more pre- 
cisely, to ip ~ 1/V12-7T ~ 1/6. Hence solution l|33(l is a good approximation 
only when the scalar field exceeds the Planckian value. This does not mean, 
however, that we require a theory of nonperturbative quantum gravity. Nonper- 
turbative quantum gravity effects become relevant only if the curvature or the 
energy density reach the Planckian values. However, even for very large values 
of the scalar field they can still remain in the sub-Planckian domain. In fact, 
considering a massive homogeneous field with negligible kinetic energy we infer 
that the energy density reaches the Planckian value for ip ~ m . Therefore, 
if m <g; 1, then for mT x > ip > 1 we can safely disregard nonperturbative 
quantum gravity effects. 
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According to l)33[) the scalar field decreases linearly with time after joining 
the attractor. During the inflationary stage 

p ~ — e + m 2 /127r. 

So when the potential energy density ~ m 2 ip 2 , which dominates the total energy 
density, drops to m 2 , inflation is over. At this time the scalar field is of order 
unity (in Planckian units). 

Let us determine the time dependence of the scale factor during inflation. 
Substituting (|33l) into (|25[) and neglecting the kinetic term, we obtain a simple 
equation which is readily integrated to yield 

a (t) ~ a f exp ( (t f - tf) ~ a, exp ( (_El±KM {t _ u) \ ^ (3g) 
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where <2j and Hi are the initial values of the scale factor and the Hubble pa- 
rameter. Note that the Hubble constant H (t) ~ y '4m [Zrrup (t) also linearly 
decreases with time. It follows from l|33[) that inflation lasts for 

At ~ t f - U ~ VUTtfai/m) . (36) 

During this time interval the scale factor increases 

^ ~ exp (27TVJ?) (37) 

times. The results obtained are in good agreement with the previous rough 
estimates (11811 and l|19ll . Inflation lasts more than 75-efolds if the initial value of 
the scalar field, tp+, is four times larger than the Planckian value. To obtain an 
estimate for the largest possible increase of the scale factor during inflation, let 
us consider a scalar field of mass 10 13 GeV. The maximal possible value of the 
scalar field for which we still remain in the sub-Planckian domain is ifi ~ 10 6 , 
and hence 

~exp(l0 12 ). (38) 

/ max 

Thus, the actual duration of the inflationary stage can massively exceed the 75 
c-folds needed. In this case our universe would constitute only a very tiny piece 
of an incredibly large homogeneous domain which originated from one causal 
region. The other important feature of inflation is that the Hubble constant 
decreases only by a factor 10~ 6 , while the scale factor grows by the tremendous 
amount given in 1)38)1 . that is, 

H; af 
Hf a t 

Graceful exit and afterwards. After the field drops below the Planckian 
value it begins to oscillate. To determine the attractor behavior in this regime 
we note that 

^ 2 + m V = ^-H 2 (39) 
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and use the Hubble parameter H and the angular variable 9, defined via 




rrup = \hj-H cos#, (40) 

V 47T 

as the new independent variables. It is convenient to replace equation (|27|l by 
a system of two first order differential equations for H and 9: 

H = -3H 2 sin 2 9, (41) 
3 

9 = -m--Hsm29, (42) 

where a dot denotes the derivative with respect to physical time t. The sec- 
ond term on the right hand side in (|42|l describes oscillations with decaying 
amplitude, as is evident from (|41|l . Therefore, neglecting this term we obtain 

9 ~ -mt + a, (43) 

where the constant phase a can be set to zero. Thus, the scalar field oscillates 
with frequency uj ~ m. After substituting 9 ~ —mt into (|41|l . we obtain a readily 
integrated equation with solution 

where a constant of integration is removed by a time shift. This solution is 
applicable only for mt ^> 1. Therefore the oscillating term is small compared 
to unity and the expression on the right hand side in (|44|l can be expanded in 
powers of (mt)~ . Substituting 143|) and (|44|l into the second equation in (|40l) . 
we obtain 

. , cos (mt) ( sm(2mt)\ /. , , 

The time dependence of the scale factor can easily be derived by integrating 

a oc t^ ( 1 - C ° S(2 ? t } ~ -Al + O Umty 3 )) . (46) 
V 6m 2 t 2 2Am 2 t 2 V ' J ) v ' 

Thus, in the leading approximation (up to decaying oscillating corrections), the 
universe expands like a matter-dominated universe with zero pressure. This 
is not surprising because an oscillating homogeneous field can be thought of 
as a condensate of massive scalar particles with zero momenta. Although the 
oscillating corrections are completely negligible in the expressions for a (t) and 
H (t) , they must nevertheless be taken into account when we calculate the 
curvature invariants. For example, the scalar curvature is 

^--^2 (l + 3cos(2mt) + o((™t) _1 )) (47) 
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(compare to R = —4/3t 2 in a matter-dominated universe). 

We have shown that inflation with a smooth graceful exit occurs naturally in 
models with classical massive scalar fields. If the mass is small compared to the 
Planck mass, the inflationary stage lasts long enough and is followed by a cold 
matter-dominated stage. This cold matter, consisting of heavy scalar particles, 
must finally be converted to radiation, baryons and leptons. We will see later 
that this can easily be achieved in a variety of ways. 



4.2 General potential: slow-roll approximation 

Equation 1(241) for a massive scalar field in an expanding universe coincides with 
the equation for a harmonic oscillator with a friction term proportional to the 
Hubble parameter H. It is well known that a large friction damps the initial ve- 
locities and enforces a slow-roll regime in which the acceleration can be neglected 
compared to the friction term. Because for a general potential H tx \fe ~ Vv, 
we expect that for large values of V the friction term can also lead to a slow-roll 
inflationary stage, where ip is negligible compared to 3H(p. Omitting the (p term 
and assuming that (p 2 <C V, equations 124(1 and (|25|l simplify to 



:»f;- + r^o. l^)-\Y v ^- (48) 



Taking into account that 



din a (ilna V^dlna 

- - - ~ v — - — - - 



dt dip 3H dip 

equations (|4*5|l give 



■^.M (49) 
dip 



and hence 



/ r<Pi v \ 
a (ip) ~ a, exp I 8tt / ji~d<p I . (50) 

This approximate solution is valid only if the slow-roll conditions 

|^ 2 |«|F|, \$\ « 3H<p ~ \V V \ , (51) 

used to simplify equations l|24|) and 1)25(1 . are satisfied. With the help of equa- 
tions 1(48(1 , they can easily be recast in terms of requirements on the derivatives 
of the potential itself: 



v 



2 



« 1, 



V 



< 1. (52) 



For a power-law potential, V = (l/n)X(p n , both conditions are satisfied for 
\(p\ ^> 1. In this case the scale factor changes as 



a (cp (t)) ~ a, exp (ip 2 - ip 2 (t))\ 



(53) 
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It is obvious that the bulk of the inflationary expansion takes place when the 
scalar field decreases by a factor of a few from its initial value. However, we are 
interested mainly in the last 50-70 e-folds of inflation because they determine the 
structure of the universe on present observable scales. The detailed picture of 
the expansion during these last 70 e-folds depends on the shape of the potential 
only within a rather narrow interval of scalar field values. 

Problem 7. Find the time dependence of the scale factor for the power law 
potential and estimate the duration of inflation. 

Problem 8. Verify that for a general potential V the system of equations 
(|24|) , 1251) can be reduced to the following first order differential equation: 

d V o ft 2\ ft , V ,a 



where 



s-^P-»n i+ #,l- < 54) 



1 47T dip 



3 V 3 d In a 

Assuming that V v /V — + as \ip\ — ► oo, draw the phase diagram and ana- 
lyze the behavior of the solutions in different asymptotic regions. Consider 
separately the case of the exponential potential. What is the physical 
meaning of the solutions in the regions corresponding to y > 1? 

After the end of inflation the scalar field begins to oscillate and the universe 
enters the stage of deceleration. Assuming that the period of oscillation is 
smaller than the cosmological time, let us determine the effective equation of 
state. Neglecting the expansion and multiplying equation l|24|) by tp, we obtain 

00)' - ip 1 + <pV <v ~ 0. (55) 

As a result of averaging over a period, the first term vanishes and hence (pi 2 ) — 
(<pV iV> ) . Thus, the averaged effective equation of state for an oscillating scalar 
field is 

P, >V^)-(2V) 

e- (pV v ) + (2V) ■ (56) 

It follows that for V cx ip n we have w ~ (n — 2) / (n + 2) . For an oscillating 
massive field (n = 2) we obtain w ~ in agreement with our previous result. 
In the case of a quartic potential (n = 4) , the oscillating field mimics an ultra- 
relativistic fluid with w ~ 1/3. 

In fact, inflation can continue even after the end of slow-roll. Considering 
the potential which behaves as 

for 1 > \ip\ ^> ip c (see Fig. 4), we infer from l|56|l that w — * — 1. This is easy 
to understand. In the case of a convex potential, an oscillating scalar field 



16 




spends most of the time near the potential walls where its kinetic energy is 
negligible and hence the main contribution to the equation of state comes from 
the potential term. 

Problem 9. Which general conditions must a potential V satisfy to provide 
a stage of fast oscillating inflation? How long can such inflation last and 
why is it not very helpful for solving the initial conditions problem? 

5 Preheating and reheating 

The theory of reheating is far from complete. Not only the details, but even 
the overall picture of inflaton decay depend crucially on the underlying particle 
physics theory beyond the Standard Model. Because there are so many possible 
extensions of the Standard Model, it does not make much sense to study the 
particulars of the reheating processes in each concrete model. Fortunately we 
are interested only in the final outcome of reheating, namely, in the possibility 
of obtaining a thermal Friedmann universe. Therefore, to illustrate the physical 
processes which could play a major role we consider only simple toy models. The 
relative importance of the different reheating mechanisms cannot be clarified 
without an underlying particle theory. However, we will show that all of them 
lead to the desired result. 

5.1 Elementary theory 

We consider an inflaton field ip of mass m coupled to a scalar field \ an d a 
spinor field ip. Their simplest interactions are described by three leg diagrams 
(Fig. 5), which correspond to the following terms in the Lagrangian: 

AL int = -gipx 2 - h<fii>i>. (57) 

We have seen that these kinds of couplings naturally arise in gauge theories 
with spontaneously broken symmetry, and they are enough for our illustrative 
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Figure 5: 



purposes. To avoid a tachyonic instability we assume that \gtp\ is smaller than 
the squared "bare" mass rnt. The decay rates of the inflaton field into XX~ 
and i^ip — pairs are determined by the coupling constants g and h respectively. 
They can easily be calculated and the corresponding results are cited in every 
book on particle physics: 



2 i 2 

r x = rfo>-> = r^ = r( v ->^) = -^. (58) 

8wm on 



Let us apply these results in order to calculate the decay rate of the inflaton. 
As we have noted, an oscillating homogeneous scalar field can be interpreted as 
a condensate of heavy particles of mass to "at rest," that is, their 3- momenta k 
are equal to zero. Keeping only the leading term in (|45|l . we have 

ip (t) ~ $ (i) cos (mi) , (59) 

where $ (t) is the slowly decaying amplitude of oscillations. The number density 
of (^-particles can be estimated as 



2 to 



i: = — = — Up 2 + mV) - T;m^ 2 . (60) 



This number is very large. For example, for m ~ 10 13 GeV, we have n v ~ 10 92 
cm" immediately after the end of inflation, when $ ~ 1 in Planckian units. 

One can show that quantum corrections do not significantly modify the 
interactions l|57(l only if g < m and h < to 1 / 2 . Therefore, for m -C mpi, the 
highest decay rate into %— particles, T x ~ m, is much larger than the highest 
possible rate for the decay into fermions, ~ to 2 . If g ~ m, then the lifetime of 
a (p— particle is about T" 1 ~ to" 1 and the inflaton decays after a few oscillations. 
Even if the coupling is not so large, the decay can still be very efficient. The 
reason is that the effective decay rate into bosons, r e g, is equal to T x , given 
in (JHSJ, only if the phase space of \— particles is not densely populated by 
previously created \— particles. Otherwise T e g can be made much larger by the 
effect of Bose condensation. This amplification of the inflaton decay is discussed 
in the next section. 

Taking into account the expansion of the universe, the equations for the 
number densities of the ip— and x~ particles can be written as 

1 d (a 3 n v ) 1 d (a 3 n x ) 

-t. = -r cff n v ; = 2T cS n v , (61) 

a A dt a 6 dt 
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where the coefficient two in the second equation arises because one (^—particle 
decays into two x~ particles. 

Problem 10. Substituting l|6U[) into the first equation in Q6ip. derive the 
approximate equation 

<p + (3H + r cff ) ip + m 2 Lp ~ 0, (62) 

which shows that the decay of the inflaton amplitude due to particle pro- 
duction may be roughly taken into account by introducing an extra fric- 
tion term T e gip. Why is this equation applicable only during the oscillatory 
phase? 

5.2 Narrow resonance 

The domain of applicability of elementary reheating theory is limited. Bose 
condensation effects become important very soon after the beginning of the in- 
flaton decay. Because the inflaton particle is "at rest," the momenta of the two 
produced x~ particles have the same magnitude k but opposite directions. If 
the corresponding states in the phase space of %— particles are already occu- 
pied, then the inflaton decay rate is enhanced by a bose factor. The inverse 
decay process XX —* <P can also take place. The rates of these processes are 
proportional to 

\{n v - l,n k + l,n_ k + 1| a^at^a^ \n v , n k , n_ k ) | 2 = (n k + 1) (n_ k + l)n v 
and 

| (n v + 1, n k - 1, n_ k - 1| aja k aZ k \n v , n k , n_ k ) | 2 = n k n_ k (n v + 1) 

respectively, where are the creation and annihilation operators for x~ parti- 
cles and n± k are their occupation numbers. To avoid confusion the reader must 
always distinguish the occupation numbers from the number densities keeping in 
mind that the occupation number refers to a density per cell of volume (2tt) (in 
the Planckian units) in the phase space, while the number density is the number 
of particles per unit volume in the three-dimensional space. Taking into account 
that n k = ?i_ k = Ufe and n v 3> 1, we infer that the number densities n v and n x 
satisfy (JSTJ, where 

r cff ~r x (l + 2n fc ). (63) 

Given a number density n x , let us calculate rife. A (p-particle "at rest" de- 
cays into two x~ particles, both having energy m/2. Because of the interaction 
term the effective squared mass of the x~ particle depends on the value 

of the inflaton field and is equal to + 2gip (t) . Therefore the corresponding 
3-momentum of the produced x~ particle is given by 

k= ((f) 2 -™x- 2 ^(*)) 12 > ( 64 ) 
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where we assume that to 2 + 2gip <C m 2 . The oscillating term, 

gip ~ g$> cos {mi) , 

leads to a "scattering" of the momenta in phase space. If g$ <C m 2 /8, then all 
particles are created within a thin shell of width 

Ak~m (^^j < m (65) 

located near the radius &o — to/2 (Fig. 6a). Therefore 

n fe - m / 2 ~ o ~ £ = -. (66) 

1 (47rfc 2 Afc) / (2ir) 3 rng<& 9 n v v ' 

The occupation numbers rife exceed unity, and hence, the Bose condensation 
effect is essential only if 

Taking into account that at the end of inflation $ ~ 1 , we infer that the occupa- 
tion numbers begin to exceed unity as soon as the inflaton converts a fraction g 
of its energy to \~ particles. The derivation above is valid only for g$> <C m 2 /8. 
Therefore, if m ~ 10~ 6 , then at most a fraction g ~ m 2 ~ 10~ 12 of the in- 
flaton energy can be transferred to x~ particles in the regime where rife < 1. 
Thus, the elementary theory of reheating, which is applicable for rife -C 1, fails 
almost immediately after the beginning of reheating. Given the result in (1661) . 
the effective decay rate becomes 

r,„-4-(i + (68) 
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where we have used equation 158fl for T x . Substituting this expression into the 
second equation in (|6 II) . we obtain 



1 d (a 3 n x ) g 2 ( ^ 27r 2 $ n 



1 + -)n tp , (69) 



a 3 dN 2m 2 \ g 

where N = mt /2tt is the number of inflaton oscillations. Let us neglect for a 
moment the expansion of the universe and disregard the decrease of the inflation 
amplitude due to particle production. In this case $ = const and for rik ^S> 1 
equation (|69|l can be easily integrated. The result is 

/7T 2 (7$ \ 

n x cx exp ( — ^-N J cx exp (2iruN) , (70) 

where fi = ■ng^j (2m 2 ) is the parameter of instability. 

Problem 11. Derive the following equation for the Fourier modes of the field 
X in Minkowski space: 

Xk + (fc 2 + m 2 + 2 ff $ cos mt) X k = 0. (71) 

Reduce it to the well-known Mathieu equation and, assuming that m 
m x 2 \g<f> | , investigate the narrow parametric resonance. Determine the 
instability bands and the corresponding instability parameters. Compare 
the width of the first instability band with (|65|l . Where is this band lo- 
cated? The minimal value of the initial amplitude of Xk is due to vacuum 
fluctuations. The increase of Xk with time can be interpreted as the pro- 
duction of x~ particles by the external classical field (p, with n x cx \xk\ 2 ■ 
Show that in the center of the first instability band, 

n x « exp (j^-Nj , (72) 

where N is the number of oscillations. Compare this result with l|70|l and 
explain why they are different by a numerical factor in the exponent. Thus, 
Bose condensation can be interpreted as a narrow parametric resonance 
in the first instability band, and vice versa. Give a physical interpretation 
of the higher order resonance bands in terms of particle production. 

Using the results of this problem we can reduce the investigation of the 
inflaton decay due to the coupling 

A£ mt - -iffW, (73) 

to the case studied above. In fact, the equation for a massless scalar field Xi 
coupled to the inflaton ip = $ cos mi, takes the form 

Xk + (fc 2 + g 2 $ 2 cos 2 mt) xk = 0, (74) 
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which coincides with equation l|71|l for to 2 — 2<?$ after substitutions <? 2 $ 2 — > 
4g<I> and to — > to/2. Thus, the two problems are mathematically equivalent. 
Using this observation and making the corresponding replacements in formula 
l|72|) . we immediately find that 

n x « exp [~j^Nj . (75) 

The condition for narrow resonance is <?$ -C to and the width of the first 
resonance band can be estimated from H65fl as Afc ~ to (g 2 $ 2 /to 2 ) . 

In summary, we have shown that even for a small coupling constant the 
elementary theory of reheating must be modified to take into account the Bose 
condensation effect, and that this can lead to an exponential increase of the 
reheating efficiency. 

Problem 12. Taking a few concrete values for g and to, compare the re- 
sults of the elementary theory with those obtained for narrow parametric 
resonance. 

So far we have neglected the expansion of the universe, the backreaction of 
the produced particles and their rescatterings. All these effects work to suppress 
the efficiency of the narrow parametric resonance. The expansion shifts the 
momenta of the previously created particles and takes them out of the resonance 
layer (Fig. 6a). Thus, the occupation numbers relevant for Bose condensation 
are actually smaller than what one would expect according to the naive estimate 
Ij66(l . If the rate of supply of newly created particles in the resonance layer is 
smaller than the rate of their escape, then rik < 1 and we can use the elementary 
theory of reheating. The other important effect is the decrease of the amplitude 
$ (t) due to both the expansion of the universe and particle production. Because 
the width of the resonance layer is proportional to it becomes more and more 
narrow. As a result the particles can escape from this layer more easily and they 
do not stimulate the subsequent production of particles. The rcscattcring of the 
X— particles also suppresses the resonance efficiency by removing particles from 
the resonance layer. Another effect is the change of the effective inflaton mass 
due to the newly produced \~ particles; this shifts the center of the resonance 
layer from its original location. 

To conclude, narrow parametric resonance is very sensitive to the interplay of 
different complicating factors. It can be fully investigated only using numerical 
methods. From our analytical consideration we can only say that the inflaton 
field probably decays not as "slowly" as in the elementary theory, but not as 
"fast" as in the case of pure narrow parametric resonance. 

5.3 Broad resonance 

So far we have considered only the case of a small coupling constants. Quantum 

1/2 

corrections to the Lagrangian are not very crucial if g < m and g < (m/$) 
They can therefore be ignored when we consider inflaton decay in the strong 
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coupling regime: m > g > m 2 /$ for the three-leg interaction and (m/^) 1 ^ 2 > 
g > to/$ for the quartic interaction J73J). In this case the condition for narrow 
resonance is not fulfilled and we cannot use the methods above. Perturbative 
methods fail because the higher order diagrams, built from the elementary dia- 
grams, give comparable contributions. Particle production can be treated only 
as a collective effect in which many inflaton particles participate simultaneously. 
We have to apply the methods of quantum field theory in an external classical 
background — as in Problem 1.11. 

Let us consider quartic interaction (|73|1 . First, we neglect the expansion of 
the universe. For ^> to the mode equation (see J71J)): 

X k + uj 2 (t) X k = 0, (76) 

where 

w(t)= (fc 2 +5 2 $ 2 cos 2 mi) 1/2 , (77) 

describes a broad parametric resonance. If the frequency lu (t) is a slowly varying 
function of time or, more precisely, \uj\ <C uj 2 , equation l|76[) can be solved in the 
quasiclassical (WKB) approximation: 

Xk oc -i= exp J Lddt^J . (78) 

In this case the number of particles, n x ~ S-xl^t i s an adiabatic invariant and 
is conserved. For most of the time the condition |w| <g; uj 2 is indeed fulfilled. 
However, every time the oscillating inflaton vanishes at tj — mT 1 (j + 1/2) 7r, 
the effective mass of the x field, proportional to |cos (mt) \ , vanishes. It is shortly 
before and after tj that the adiabatic condition is strongly violated: 

\uj\ _ mg 2 <5> 2 |cos (mi) sin (mt)\ > ^ 
u 2 (fc 2 + g 2 $ 2 cos 2 (mt)f /2 

Considering a small time interval At m~ 1 in the vicinity of tj, we can rewrite 
this condition as 

^Atl + (At/AQ 2 ) 

where 

At, ~ (g<$> m y 1/2 = — (g<S>/ m y 1/2 . (81) 

TO 

It follows that the adiabatic condition is broken only within short time intervals 
At ~ At, near tj and only for modes with 

k < h ~ At' 1 ~m(g<i>/m) 1/2 . (82) 

Therefore, we expect that particles with the corresponding momenta are 
created only during these time intervals. It is worth noting that the momen- 
tum of the created particle can be larger than the inflaton mass by the ratio 
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{g^/m 2 ) 1 ' 2 > 1; the X~ particles are produced as a result of a collective process 
involving many inflaton particles. This is the reason why we cannot describe 
the broad resonance regime using the usual methods of perturbation theory. 

To calculate the number of particles produced in a single inflaton oscillation 
we consider a short time interval in the vicinity of tj and approximate the cosine 
in (|76[) by a linear function. Equation <|76l) then takes the form 

^ +(k 2 +t 2 ) Xk = 0, (83) 



dr 2 



where the dimensionless wavenumbcr k = k/k* and time r = (t — tj) /At* 
have been introduced. In terms of the new variables the adiabaticity condition 
is broken at |r| < 1 and only for k < 1. It is remarkable that the coupling 
constant g, the mass and the amplitude of the inflaton enter explicitly only in 
k 2 . The adiabaticity violation is largest for k = 0. In this case the parameters 
g, $ and m drop from equation Q83|l and the amplitude x«=o changes only by 
a numerical, parameter-independent factor as a result of passing through the 
nonadiabatic region at |r| < 1. Because the particle density n is proportional to 
\x\ , its growth from one oscillation to the next can written as 

= exp (27T/x fe=0 ) , (84) 



where the instability parameter fik=o does not depend on g, $ and m. For 
modes with k ^ 0, the parameter /ik^o is a function of k = k/k*. In this case 
the adiabaticity is not violated as strongly as for the k = mode, and hence 
/ifc^o is smaller than fj,k=o- To calculate the instability parameters we have to 
determine the change of the amplitude \ i n passing from the r < — 1 region to 
the t > 1 region. This can be done using two independent WKB solutions of 
equation (|83|l in the asymptotic regions |r| 1: 

X± = (K 2 + 1 T 2 ) i/4 exp { ±l J ^+^ dT ) - kP^exp (±!0 . (85) 

After passing through the nonadiabatic region the mode A + x+ becomes a mix- 
ture of the modes x+ an d X- i th&t is, 

A +x+ ^ B + x+ + C + x-, (86) 

where A + ,B + and C+ are the complex constant coefficients. Similarly, for the 
mode A^x~i we have 

A„x- -^B_x-+C-x+- (87) 

Drawing an analogy with the scattering problem for the inverse parabolic po- 
tential, we note that the mixture arises due to an overbarrier reflection of the 
wave. The reflection is most efficient for the waves with k = which "touch" 
the top of the barrier. 



24 



The quasiclassical solution is valid in the complex plane for |t| ^> 1. Travers- 
ing the appropriate contour t = |t| e lv in the complex plane from t <C — 1 to 
t > 1, we infer from (JSSJ), (HSU) and JSZJl that 

B ± = T ie"7 K2 yl ± . (88) 

The coefhcients C± are not determined in this method. To find them we use 
the Wronskian 

W = XX*~XX\ (89) 

where x is an arbitrary complex solution of equation l|83|) . Taking the derivative 
of W and using (|83[l to express x m terms of x, we find 

W = (90) 

and hence W = const. From this we infer that the coefficients A 7 B and C in 
1)86(1 and H87|l satisfy the "probability conservation" condition 

|C ± | 2 -|5 ± | 2 = L4 ± | 2 . (91) 
Substituting B from l|88|l . we obtain 



C± = VT+e 3 ™ 1 \A± | e M± , (92) 

where the phases a± remain undetermined. 

At |t| ^ 1 the modes of field x satisfy the harmonic oscillator equation with 
a slowly changing frequency u> oc |r| . In quantum field theory the occupation 
number nk in the expression for the energy of the harmonic oscillator, 

e k =uj(n k + 1/2) , (93) 

is interpreted as the number of particles in the corresponding mode k. In the 
adiabatic regime (|t| ^> 1) this number is conserved and it changes only when 
the adiabatic condition is violated. Let us consider an arbitrary initial mixture 
of the modes x+ an d X-- After passing through the nonadiabatic region at 
t ~ tj, it changes as 

X 3 - A +x+ + A^x- - X 3+1 = (B+ + CL) x+ + (B- + C+) X - (94) 
Taking into account that 

n+i = -^|x| 2 , (95) 

we see that as a result of this passage the number of particles in the mode k 
increases 



1/2 \ ^|x J+1 | \B+ + C-\ 2 + \B-+C+\ 2 



nJ + l/2 J k w | x i| 2 " \A + \ 2 + \AJ 2 



(96) 
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times, where we have averaged \x\ over the time interval m 1 > t > u> 1 . With 
B and C from (|88|l and i|92|) , this expression becomes 

V nP + 1/2 ) k V J \A+\ 2 + \A_\ 2 

Problem 13. Verify the result l|97|l and explain the origin of the phase 8. 

In the vacuum initial state n,}. = but the amplitude of the field \ does not 
vanish because of the existence of vacuum fluctuations; we have |^4+| ^ and 
|j4°_| 2 = 0. It follows from the "probability conservation" condition that 

L4 + | 2 -|A_| 2 = |^| 2 (98) 

at every moment of time. This means that as a result of particle production the 
coefficients \A + \ 2 and |A_| 2 grow by the same amount. When \A + \ becomes 
much larger than \A°^\ we have \A + \ ~ \A_ \ . Taking this into account and 
beginning in the vacuum state, we find from (|97|l that after N ^ 1 inflaton 
oscillations the particle number in mode k is 

n k ~ 1 exp (27rA**JV) , (99) 

where the instability parameter is given by 

H k ~ — In f 1 + 2e- 7 " t2 +2cos(9e^ K Vl + e- ,rK2 ) . (100) 

This parameter takes its maximal value 

/ir x = 7r _1 ln(l + V2) -0.28 

for k — and 8 = 0. In the interval — tt < 8 < 7r we find that /ik=o is positive if 
37r/4 > 9 > — 37r/4 and negative otherwise. Thus, assuming random 8, we con- 
clude that the particle number in every mode changes stochastically. However, 
if all 8 are equally probable, then the number of particles increases three quar- 
ters of the time and therefore it also increases on average, in agreement with 
cntropic arguments. The net instability parameter, characterizing the average 
growth in particle number, is obtained by skipping the cos#-term in <|100|) : 

fi k ~ i.m(l + 2 e -™ 2 ) . (101) 

With slight modifications the results above can be applied to an expanding 
universe. First of all we note that the expansion randomizes the phases 8 and 
hence the effective instability parameter is given by l|101|) . For particles with 
physical momenta k < ksg/s/n, the instability parameter p, k can be roughly 
estimated by its value at the center of the instability region, p, k =o = ln3/27r ~ 
0.175. To understand how the expansion can influence the efficiency of broad 
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resonance, it is again helpful to use the phase space picture. The particles 
created in the broad resonance regime occupy the entire sphere of radius A;*/ ' ypH 
in phase space (see Fig. 6b). During the passage through the nonadiabatic 
region the number of particles in every cell of the sphere, and hence the total 
number density, increases on average exp (2ir ■ 0.175) ~ 3 times. At the stage 
when inflaton energy is still dominant, the physical momentum of the created 
particle decreases in inverse proportion to the scale factor (fc oc a -1 ), while 
the radius of the sphere shrinks more slowly, namely, as < I> 1 / 2 oc < -1 / 2 oc a -3 / 4 . 
As a result, the created particles move away from the boundary of the sphere 
towards its center where they participate in the next "act of creation," enhancing 
the probability by a bose factor. Furthermore, expansion also makes broad 
resonance less sensitive to rescattering and backreaction effects. These two 
effects influence the resonance efficiency by removing those particles which are 
located near the boundary of the resonance sphere. Because expansion moves 
particles away from this region, the impact of these effects is diminished. Thus, 
in contrast to the narrow resonance case, expansion stabilizes broad resonance 
and at the beginning of reheating it can be realized in its pure form. 

Taking into account that the initial volume of the resonance sphere is about 

kl ~ to 3 (g<P /mf /2 , 

we obtain the following estimate for the ratio of the particle number densities 
after TV inflaton oscillations: 



% _ fc 3 exp (2irn k=0 N) ^ ^1/2-3/2 ., 
n M m<F 2 



o 



where $o ~ O (1) is the value of the inflaton amplitude after the end of inflation. 
Since in the adiabatic regime the effective mass of the x~ particles is of order 
<7$, where $ decreases in inverse proportion to N, we also obtain an estimate 
for the ratio of the energy densities: 

£* ^ UhlhL „ m -i/2 5 6/2j V -i 3 Jv_ (103) 

The formulae above fail when the energy density of the created particles begins 
to exceed the energy density stored in the inflaton field. In fact, at this time, 
the amplitude $ (t) begins to decrease very quickly because of the very efficient 
energy transfer from the inflaton to the \~ particles. Broad resonance is cer- 
tainly over when $ (t) drops to the value $ r ~ m/g, and we enter the narrow 
resonance regime. For the coupling constant to 1 / 2 > g > 0(1) m, the num- 
ber of the inflaton oscillation N r in the broad resonance regime can be roughly 
estimated using the condition e x <~ e v : 

N r - (0.75 to 2) log 3 m" 1 . (104) 

As an example, if m ~ 10 13 GcV, we have N r ~ 10 to 25 for a wide range of the 
coupling constants 10~ 3 > g > 10~ 6 . Taking into account that the total energy 
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decays as m 2 (<&o/A0 , we obtain 





) 



2 



e. 



(105) 



e y + e 



2 {^/N r f 



that is, the energy still stored in the inflaton field at the end of broad resonance 
is only a small fraction of the total energy. In particular, for m ~ 10 13 GeV, 
this ratio varies in the range 1CP 6 to O (1) depending on the coupling constant 



Problem 14. Investigate inflaton decay due to the three-leg interaction g 2 ipx 2 
in the strong coupling regime: m > g > m 2 /$. 

5.4 Implications 

It follows from the above considerations that broad parametric resonance can 
play a very important role in the preheating phase. During only 15 to 25 
oscillations of the inflaton, it can convert most of the inflaton energy into other 
scalar particles. The most interesting aspect of this process is that the effective 
mass and the momenta of the particles produced can exceed the inflaton mass. 
For example, for m ~ 10 14 GeV, the effective mass = g&\cos(mt)\ can 
be as large as 10 16 GeV. Therefore, if the \~ particles are coupled to bosonic 
and fermionic fields heavier than the inflaton, then the inflaton may indirectly 
decay into these heavy particles. This brings Grand Unification scales back into 
play. For instance, even if inflation ends at low energy scales, preheating may 
rescue the GUT baryogenesis models. Another potential outcome of the above 
mechanism is the far-from-equilibrium production of topological defects after 
inflation. Obviously their numbers must not conflict with observations and this 
leads to cosmological bounds on admissible theories. 

If, after the period of broad resonance, the slightest amount of the infla- 
ton remained — given by 1)105(1 — it would be a cosmological disaster. Since 
the inflaton particles are nonrelativistic, if they were present in any substantial 
amount, they would soon dominate and leave us with a cold universe. Fortu- 
nately, these particles should easily decay in the subsequent narrow resonance 
regime or as a result of elementary particle decay. These decay channels thus 
become necessary ingredients of the reheating theory. 

The considerations of this section do not constitute a complete theory of 
reheating. We have studied only elementary processes which could play a role 
in producing a hot Friedmann universe. The final outcome of reheating must 
be matter in thermal equilibrium. The particles which are produced in the 
preheating processes are initially in a highly non-equilibrium state. Numerical 
calculations show that as a result of their scatterings they quickly reach local 
thermal equilibrium. Parametrizing the total preheating and reheating time in 
terms of the inflaton oscillations number Nt, we obtain the following estimate 
for the reheating temperature: 
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where N is the effective number degrees of freedom of the light fields at T <~ Tr. 
Assuming that N T ~ 10 6 , and taking N <~ 10 2 and m ~ 10 13 GeV, we obtain 
Tr ~ 10 12 GeV. This does not mean, however, that we can ignore physics 
beyond this scale. As we have already pointed out, non-equilibrium preheating 
processes can play a nontrivial role. 

Reheating is an important ingredient of inflationary cosmology. We have 
seen that there is no general obstacle to arranging successful reheating. A par- 
ticle theory should be tested on its ability to realize reheating in combination 
with baryogenesis. In this way, cosmology enables us to preselect realistic par- 
ticle physics theories beyond the Standard Model. 

6 "Menu" of scenarios 

All we need for successful inflation is a scalar condensate satisfying the slow-roll 
conditions. Building concrete scenarios then becomes a "technical" problem. 
Involving two or more scalar condensates, and assuming them to be equally rel- 
evant during inflation, extends the number of possibilities, but simultaneously 
diminishes the predictive power of inflation. This especially concerns cosmolog- 
ical perturbations, which are among the most important robust predictions of 
inflation. Because inflation can be falsified experimentally (or more accurately, 
observationally) only if it makes such predictions, we consider only simple sce- 
narios with a single inflaton component. Fortunately all of them lead to very 
similar predictions which differ only slightly in the details. This makes the signif- 
icance of a unique scenario, the one actually realized in nature, less important. 
The situation here is very different from particle physics, where the concrete 
models are as important as the ideas behind them. This does not mean we do 
not need the correct scenario; if one day it becomes available, we will be able 
to verify more delicate predictions of inflation. However, even in the absence of 
the true scenario, we can nonetheless verify observationally the most important 
predictions of the stage of cosmic acceleration. The purpose of this section is to 
give the reader a very brief guide to the "menu of scenarios" discussed in the 
literature. 

Inflaton candidates. The first question which naturally arises is "what is 
the most realistic candidate for the inflaton field?" . There are many because the 
only requirement is that this candidate imitates a scalar condensate in the slow- 
roll regime. This can be achieved by a fundamental scalar field or by a fcrmionic 
condensate described in terms of an effective scalar field. This, however, does 
not exhaust all possibilities. The scalar condensate can also be imitated entirely 
within the theory of gravity itself. Einstein gravity is only a low curvature limit 
of some more complicated theory whose action contains higher powers of the 
curvature invariants, for example, 

S = ( (R + aR 2 + f3R^R^ + 1 R 3 + ...) ^jdx. (107) 

l07T J 

The quadratic and higher order terms can be cither of fundamental origin or they 
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can arise as a result of vacuum polarization. The corresponding dimensional 
coefficients in front of these terms are likely of Planckian size. The theory 
with action (|107fl can provide us with inflation. This can easily be understood. 
Einstein gravity is the only metric theory in four dimensions where the equations 
of motion are second order. Any modification of the Einstein action introduces 
higher derivative terms. This means that, in addition to the gravitational waves, 
the gravitational field has extra degrees of freedom including, generically, a spin 
zero field. 

Problem 15. Consider a gravity theory with metric g M „ and action 

5= 77T / f(R)V=9dx, (108) 

167T J 

where / (i?) is an arbitrary function of the scalar curvature R. Derive the 
following equations of motion: 

£*-5«' + (fQ>-(i9> (109) 

Verify that under the conformal transformation — > g^ v = Fg^ v , the 
Ricci tensor and the scalar curvature transform as 

R£^R£ = F- 1 ^ - F- 2 F:» - ^F- 2 F;«5Z + -F^F^F*, (110) 

-> R = F~ l R - 3F~ 2 F£ + -F~ a F ta F< a . (Ill) 
Introduce the "scalar field" 



^Vii lnF(i?)i (112) 

and show that the equations 

k - \m = mts (v>) > (us) 

coincide with (|109J) if we set F — df jdR and take the following potential 
for the scalar field: 



1 f-ROf/d R 
16tt (df/OR) 2 



y^)=ih j ,::::z ■ ( ii4 ) 



Problem 16. Study the inflationary solutions in R 2 — gravity: 

1 f /„ 1 



S = ~Wj [ R -6AP R2 )^ dX - (U5) 
What is the physical meaning of the constant Ml 
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Thus, the higher derivative gravity theory is conformally equivalent to Ein- 
stein gravity with an extra scalar field. If the scalar field potential satisfies 
the slow-roll conditions, then we have an inflationary solution in the conformal 
frame for the metric g^ v . One should not confuse, however, the conformal metric 
with the original physical metric. They generally describe manifolds with differ- 
ent geometries and the final results must be interpreted in terms of the original 
metric. In our case the use of the conformal transformation is a mathemati- 
cal tool which simply allows us to reduce the problem to one we have studied 
before. The conformal metric is related to the physical metric g^ v by a factor 
F, which depends on the curvature invariants; it does not change significantly 
during inflation. Therefore, we also have an inflationary solution in the original 
physical frame. 

So far we have been considering inflationary solutions due to the potential of 
the scalar field. However, inflation can be realized even without a potential term. 
It can occur in Born-Infeld type theories, where the action depends nonlinearly 
on the kinetic energy of the scalar field. These theories do not have higher 
derivative terms, but they have some other peculiar properties. 

Problem 17. Consider a scalar field with action 

S = J p(X,<p)^dx, (116) 

where p is an arbitrary function of ip and X = ^ (d^ipdtp^) . Verify that 
the energy-momentum tensor for this field can be written in the form 

T$ = {e + p)vPu v -p&Z, (117) 

where the Lagrangian p plays the role of the effective pressure and 

(118) 

If the Lagrangian p satisfies the condition Xdp/dX <C p for some range of 
X and ip, then the equation of state is p — e and we have an inflationary 
solution. Why is inflation not satisfactory if p depends only on XI Con- 
sider a general p (X, <p) without an explicit potential term, that is, p — > 
when X — > 0. Formulate the conditions which this function must satisfy 
to provide us with a slow-roll inflationary stage and a graceful exit. The 
inflationary scenario based on the nontrivial dependence of the Lagrangian 
on the kinetic term is called fc— inflation. 



Scenarios. The simplest inflationary scenarios can be subdivided into three 
classes. They correspond to the usual scalar field with a potential, higher deriva- 
tive gravity, and fc— inflation. The cosmological consequences of scenarios from 
the different classes are almost indistinguishable — they can exactly imitate 
each other. Within each class, however, we can try to make further distinctions 
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by addressing the questions: a) what was before inflation? and b) how does a 
graceful exit to a Friedmann stage occur? For our purpose it will be sufficient 
to consider only the simplest case of a scalar field with canonical kinetic energy. 
The potential can have different shapes, as shown in Fig. 7. The three cases pre- 
sented correspond to the so-called old, new and chaotic inflationary scenarios. 
The first two names refer to their historical origins (see "Bibliography" ) . 

Old inflation (see Fig. 7a) assumes that the scalar field arrives at the local 
minimum of the potential at ip — as a result of a supercooling of the initially 
hot universe. After that the universe undergoes a stage of accelerated expansion 
with subsequent graceful exit via bubble nucleation. It was clear from the very 
beginning that this scenario could not provide a successful graceful exit because 
all the energy released in a bubble is concentrated in its wall and the bubbles 
have no chance to collide. This difficulty was avoided in the new inflationary 
scenario, a scenario similar to a successful model in higher derivative gravity 
which had previously been invented. 

New inflation is based on a Coleman- Weinberg type potential (Fig. 7b). 
Because the potential is very fiat and has a maximum at ip = 0, the scalar 
field escapes from the maximum not via tunneling, but due to the quantum 
fluctuations. It then slowly rolls towards the global minimum where the energy 
is released homogeneously in the whole space. Originally the pre-inflationary 
state of the universe was taken to be thermal so that the symmetry was restored 
due to thermal corrections. This was a justification for the initial conditions 
of the scalar field. Later on it was realized that the thermal initial state of 
the universe is quite unlikely, and so now the original motivation for the initial 
conditions in the new inflationary model seems to be false. Instead, the universe 
might be in a "self-reproducing" regime (for more details see section 8.5). 

Chaotic inflation gives its name to the broadest possible class of potentials 
satisfying the slow roll-conditions (Fig. 7c). We have considered it in detail in 
the previous sections. The name chaotic is related to the possibility of having 
almost arbitrary initial conditions for the scalar field. To be precise, this field 
must initially be larger than the Planckian value but it is otherwise arbitrary. 
Indeed, it could have varied from one spatial region to another and, as a result, 
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the universe would have a very complicated global structure. It could be very 
inhomogencous on scales much larger than the present horizon and extremely 
homogeneous on "small" scales corresponding to the observable domain. We 
will see later that in the case of chaotic inflation, quantum fluctuations lead to 
a self-reproducing universe. 

Since chaotic inflation encompasses so many potentials, one might think it 
worthwhile to consider special cases, for example, an exponential potential. For 
an exponential potential, if the slow-roll conditions are satisfied once, they are 
always satisfied. Therefore, it describes (power-law) inflation without a graceful 
exit. To arrange a graceful exit we have to "damage" the potential. For two or 
more scalar fields the number of options increases. Thus it is not helpful here 
to go into the details of the different models. 

In the absence of the underlying fundamental particle theory, one is free to 
play with the potentials and invent more new scenarios. In this sense the situ- 
ation has changed since the time the importance of inflation was first realized. 
In fact, in the 80s many people considered inflation a useful application of the 
Grand Unified Theory that was believed to be known. Besides solving the initial 
conditions problem, inflation also explained why we do not have an overabun- 
dance of the monopoles that are an inevitable consequence of a GUT. Either 
inflation ejects all previously created monopoles, leaving less than one monopole 
per present horizon volume, or the monopoles are never produced. The same 
argument applies to the heavy stable particles that could be overproduced in the 
state of thermal equilibrium at high temperatures. Many authors consider the 
solution of the monopole and heavy particle problems to be as important as a 
solution of the initial conditions problem. We would like to point out, however, 
that the initial conditions problem is posed to us by nature, while the other 
problems are, at present, not more than internal problems of theories beyond 
the Standard Model. By solving these extra problems, inflation opens the door 
to theories that would otherwise be prohibited by cosmology. Depending on 
one's attitude, this is either a useful or damning achievement of inflation. 

De Sitter solution and inflation. The last point we would like to make 
concerns the role of a cosmological constant and a pure de Sitter solution for 
inflation. We have already said that the pure de Sitter solution cannot provide 
us with a model with a graceful exit. Even the notion of expansion is not 
unambiguously defined in de Sitter space. We saw in section 1.3.6 that this space 
has the same symmetry group as Minkowski space. It is spatially homogeneous 
and time translation invariant. Therefore any space-like surface is a hypersurface 
of constant energy. To characterize an expansion we can use not only k = 
0, ±1 Friedmann coordinates but also, for example, "static coordinates" (see 
Problem 2.7), which describe an expanding space outside the event horizon. 
In all these cases the three-geometries of constant time hypersurfaces are very 
different. These differences, however, simply characterize the different slicings 
of the perfectly symmetrical space and there is no obvious preferable choice for 
the coordinates. 

It is important therefore that inflation is never realized by a pure de Sitter 
solution. There must be deviations from the vacuum equation of state, which 
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Figure 8: 



finally determine the "hypersurface" of transition to the hot universe. The 
de Sitter universe is still, however, a very useful zeroth order approximation for 
nearly all inflationary models. In fact, the effective equation of state must satisfy 
the condition e + 3p < for at least 75 e-folds. This is generally possible only if 
during most of the time we have p w — e to a rather high accuracy. Therefore, 
one can use the language of constant time hypersurfaces defined in various 
coordinate systems in de Sitter space. Our earlier considerations show that the 
transition from inflation to the Friedmann universe occurs along a hypersurface 
of constant time in the expanding isotropic coordinates (n = const), but not 
along r = const hypersurface in the "static coordinates." The next question is, 
out of the three possible isotropic coordinate systems (k — 0, ±1), which must be 
used to match the de Sitter space to the Friedmann universe? Depending on the 
answer to this question, we obtain flat, open or closed Friedmann universes. It 
turns out, however, that this answer seems not to be relevant for the observable 
domain of the universe. In fact, if inflation lasts more than 75 e-folds, the 
observable part of the universe corresponds only to a tiny piece of the matched 
global conformal diagrams for de Sitter and Friedmann universes. This piece is 
located near the upper border of the conformal diagram for de Sitter space and 
the lower border for the flat, open or closed Friedmann universes (Fig. 8), where 
the difference between the hypersurfaces of constant time for flat, open or closed 
cases is negligibly small. After a graceful exit we obtain a very large domain of 
the Friedmann universe with incredibly small flatness and this domain covers 
all present observable scales. The global structure of the universe on scales 
much larger than the present horizon is not relevant for an observer — at least 
not for the next hundred billion years. In the next section we will see that the 
issue of the global structure is complicated by quantum fluctuations. These 
fluctuations are amplified during inflation and as a result the hypersurface of 
transition has "wrinkles." The wrinkles are rather small on scales corresponding 
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to the observable universe but they become huge on the very large scales. Hence, 
globally the universe is very different from the Friedmann space and the question 
about the spatial curvature of the whole universe no longer makes sense. It also 
follows that the global properties of an exact de Sitter solution have no relevance 
for the real physical universe. 
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